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We calculate fluctuation corrections to the longitudinal conductivity of disordered superconduc- 
tors subject to an external magnetic field. We derive analytic expressions that are valid in the entire 
metallic part of the temperature-magnetic field phase diagram as long as the effect of the magnetic 
field on the spin degrees of freedom of the electrons may be neglected. Our calculations are based 
on a kinetic equation approach. For the special case of superconducting films and wires in par- 
allel magnetic fields we perform a detailed comparison with results that were previously obtained 
with diagrammatic perturbation theory in the imaginary time formalism. As an application, we 
study the fluctuation conductivity of films in tilted magnetic fields with a special focus on the low- 
temperature regime. We present a detailed discussion of the phenomenon of the non-monotonous 
magnetoresistance and find that it displays a pronounced dependence on the tilting angle. 
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I. INTRODUCTION 

The theory of superconducting fluctuations has been 

the subject of intense study for many years. -'^ In the 
metallic part of the phase diagram, outside the super- 
conducting regime, Cooper pairs may form for a fi- 
nite time. The presence of these fluctuating Cooper 
pairs aSects both thermodynamic and transport prop- 
erties of the metal. The phase transition between the 
metalhc and the superconducting state may be tuned by 
temperature or by so-caUed pair-breaking mechanisms 
which lead to a partial or even complete destruction 
of superconductivity.^ Examples of pair-breakers include 
magnetic impurities in s-wave superconductors, exter- 
nal magnetic fields or a fiux penetrating a superconduc- 
tor with doubly connected geometry. Fluctuation effects 
are particularly strong for low-dimensional superconduc- 
tors and further enhanced by the presence of impurities.^ 
Detailed experimental^"^" and theoretical studies^^"^^ of 
fluctuation phenomena in superconductors have become 
available in recent years. 

The subject of this paper is the calculation of fluctu- 
ation corrections to conductivity in the metallic phase 
of disordered superconductors. The origin of this fleld 
dates back to the work of Azlamazov and Larkin.^^ These 
authors studied the direct contribution of fluctuating 
Cooper pairs to conductivity close to the transition tem- 
perature Tco, the so-called paraconductivity. Shortly af- 
terwards, additional contribiitions were discovered. 
This development went hand in hand with the study of 
different classes of diagrams in many-body perturbation 
theory. It became customary to divide the set of most rel- 
evant diagrams into three classes, the Aslamazov-Larkin 
diagram and the density of states and Maki-Thompson 
diagrams.^ 

Initially, studies focused around the vicinity of Tgo for 

vanishing or small magnetic fields. In Ref. 11, the fiuc- 
tuation conductivity was calculated for disordered su- 



perconducting films in perpendicular magnetic fields in 
the vicinity of the critical magnetic field Bc2- It was 
established that at very low temperatures superconduct- 
ing fiuctuations lead to a non-monotonous magnetoresis- 
tance (NM); close to Bci, the resistance curve displays 
a maximum as a function of the magnetic field. In an- 
other theoretical study, Ref. 12, it was found that the 
NM also exists in the vicinity of certain other classes of 
pair-breaking transitions such as for films and wires in 
parallel magnetic fields. 

Recently, a novel scheme for deriving fluctuation con- 
ductivity was introduced,^'' which is based on the Usadel 
equation. The calculation is performed in the Keldysh 
formalism to circumvent the analytic continuation nec- 
essary in the Kubo technique. For temperatures close 
to Tco and in the absence of a magnetic fleld, the Usadel 
equation has been used for the calculation of fluctuation 
conductivity before. In Ref. 20, in turn, general expres- 
sions for the fluctuation conductivity in disordered super- 
conducting films with perpendicular magnetic field were 
derived for the whole normal part of the temperature- 
magnetic field phase diagram (outside the strong fluctua- 
tion regime close to the transition line) . In this approach, 
it was possible to identify three distinct contributions to 
conductivity at the very early stages of the calculation. 
The flrst one, termed density of states correction ((5(Tdo.s)- 
is seen to be directly related to the change in the quasi- 
particle density of states. The second contribution is the 
anomalous Maki-Thompson correction {5<Jan), which is 
known from diagrammatic perturbation theory and de- 
scribes a coherent rescattering in the Cooper channel. 
The third term may be interpreted as the direct contri- 
bution of Cooper pairs to the current, and was therefore 
named supercurrent correction {5a sc)- It should be noted 
that the density of states and supercurrent contributions 
in the Usadel equation approach are in general not iden- 
tical to the contributions of the density of states and 
Aslamazov-Larkin diagrams in the conventional classifl- 
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cation. 

Let us briefly recall the origin of the NM for the per- 
pendicular magnetic field case using the language intro- 
duced in Ref. 20. We discuss the low temperature regime 
t = T/Tco <C 1 close to the (temperature-dependent) crit- 
ical field Bc2(T), so that h = {B- B,2{T))/ Bc2{T) < 1. 
For t ^ h, all corrections. Sa^og, 6(Ta„, said Sage, con- 
tribute, and the total correction to conductivity is pos- 
itive. As one moves further away from the transition 
line at fixed temperature, in the limit t ^ h, the anoma- 
lous Maki-Thompson correction becomes ineffective. The 
density of states correction to conductivity, which is neg- 
ative, and the supercurrent currection, which is positive, 
are of similar magnitude. The density of states correction 
dominates, however, leading to a net negative correction 
to conductivity. For large magnetic fields, the negative 
density of states correction is still dominant, but eventu- 
ally diminishes. The result is a non-monotonoTis magne- 
toresistancc. Close to 3^2, the results of Ref. 20 coincide 
with those obtained by Galitski and Larkin,^^ who specif- 
ically focused on this regime and used the conventional 
diagrammatic method for the calculation. It should be 
noted that the low-temperature regime is quite different 
from the well-studied case of small magnetic fields for 
T « Tco- In the latter regime. Sage is much larger in mag- 
nitude than Sa^og. The main difference is that close to 
Bc,2 Landau level quantization of the Cooper pair prop- 
agator becomes crucial. This is why the supercurrent 
correction becomes less singular and Sage and Saoos are 
of a comparable magnitude. 

In this paper, we use the Usadel equation approach to 
derive general expressions for the fluctuation conductiv- 
ity in superconductors subject to a magnetic field. We 
assume that the sample geometry is translationally in- 
variant along the direction of the electric field, while the 
sample may be confined in the transverse direction(s). 
The derived formulas are in particular applicable for su- 
perconducting wires, superconducting films in magnetic 
fields of arbitrary orientation and for cylinders threaded 
by a magnetic flux. As a speciflc application, we study in 
detail the phenomenon of the NM for films in tilted mag- 
netic fields. We focus on the low temperature regime, and 
describe the evolution of the NM as a function of the tilt- 
ing angle. While the phenomenon persists for any angle, 
there are two distinct regions, one comprising the parallel 
magnetic field case and the other one the perpendicular 
magnetic field case, for which the physical origin of the 
phenomenon as well as the magnitude of the resulting 
resistance-maximum are quite distinct. The cross-over 
between the two regimes occurs for almost parallel mag- 
netic fields. 

The theory developed in this paper is applicable in the 
limit of weak disorder, epT » 1, where ep is the Fermi 

energy and r the transport scattering time. From the 
experimental perspective, detailed low-temperature resis- 
tance measurements have been performed on weakly dis- 
ordered films in perpendicular magnetic fields, see, e.g., 
Refs. 28, 29 and 5. Measurements on films in parallel'^^''^^ 



and tilted "^^ magnetic fields exist, but focused on more 
strongly disordered films in the context of the so-called 
superconductor-insulator transition. 

For the case of superconductors in parallel magnetic 
fields, we perform a detailed comparison of our results to 
those obtained in Ref. 12 with the help of the traditional 
diagrammatic technique. We show that there is a one-to 
one correspondence between the results obtained in the 
two formalisms (up to details of the ultraviolet regular- 
ization). The mapping is not simple, however. The three 
distinct contributions identified in the kinetic equation 
approach correspond to a mixture of terms originating 
from different diagrams. This comparison is motivated 
by a discrepancy between recent results reported for films 
in perpendicular magnetic fields. Glatz, Varlamov and 
Vinokur^*'"'^'' used the traditional approach in the imag- 
inary time formalism for the calculation of the fiuctua- 
tion conductivity. The results of this study disagree with 
a number of previously reported results, including those 
for B « Bc2, Ref. 11, and the high temperature regime 
T > Tco at B = 0, Ref. 33. The technically very different 
Usadel equation approach of Ref. 20, however, confirmed 
these earlier results. The comparison performed here for 
the parallel magnetic field case demonstrates an agree- 
ment between the imaginary time formalism as worked 
out in Ref. 12 and the Usadel equation approach on the 
level of general formulas valid in the entire normal part of 
the phase diagram. This includes, in particular, the high 
temperature regime T Tco for B = 0, which lies within 
the range of applicability of all the mentioned works. 

The paper is organized as follows. In Sec. II we present 
the main results of our study. Specifically, in Sec. II A, we 
display the general formulas for the fluctuation conduc- 
tivity in superconductors subject to a magnetic fleld and 
discuss their range of applicability. In Sec. II B, we dis- 
cuss the fluctuation conductivity of a thin amorphous su- 
perconducting film in a tilted magnetic field. The exam- 
ple of the film in a tilted magnetic field constitutes a spe- 
cial application of the general formalism outlined in this 
paper. The remainder of the paper is devoted to the tech- 
nical details of the approach as well as to a comparison 
with the traditional diagrammatic technique. In Sec. Ill 
we introduce the Usadel equation approach underlying 
the calculation of fluctuation conductivity as well as the 
derivation of the results presented in Sec. II. The formal- 
ism we use is a generalization of the approach introduced 
in Ref. 20 so as to include pair-breaking effects. We out- 
line the main steps of the derivation in a condensed form 
in order to make the paper self-contained. In Sec. IV 
we specialize to the parallel magnetic field case and com- 
pare our results to those obtained in Ref. 12 using the 
traditional diagrammatic approach. The results of the 
comparison are summarized in table I. Sec. V is devoted 
to films in a tilted magnetic field. Here, we present the 
derivation of the results presented in Sec. II B. Eventu- 
ally, in Sec. VI we conclude. Some technical details of 
the calculation are relegated to two appendices. 
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II. RESULTS 



2. General expressions for fluctuation corrections 



In this section, we present the main results of our 

study. Wc split the discussion into two parts. In the first 
part, we present the results for the general theory of fluc- 
tuation transport in superconductors subject to a mag- 
netic field. In the second part, we discuss the fluctuation 
conductivity of a thin disordered film in a tilted magnetic 
field. The presentation is intended to be self-contained. 
Details of the derivation are described in Sees. Ill and V. 



A. General Results: Disordered superconductors 
in a magnetic field 

The main result of this paper are expressions for the 
fluctuation conductivity in disordered superconductors 
subject to a magnetic field. In Sec. II B below, we dis- 
cuss the case of a thin film in a tilted magnetic field as an 
application. The general results, however, are applicable 
not only to thin films, but also to several other geome- 
tries, for example wires, cylindric tubes, and nanorib- 
bons. The difference between these examples lies in the 
spectrum of superconducting fluctuations. In the follow- 
ing, we present the equation determining the fluctuation 
spectrum, Eq. (1), and write down the general results for 
the corrections to conductivity, Eqs. (3) to (9). Then, we 
briefly discuss the range of applicability. 



1. Fluctuation spectrum 

The fluctuation propagator of the superconducting or- 
der parameter field, cf. Eqs. (10) and (11) below, is di- 
agonal in the basis of eigenfunctions determined by the 
following eigenvalue equation 



D 



[V - 2ieA{T)fct>n (r) = a„(/)„ (r) . 



(1) 



This equation is similar to the single-particle Schrodinger 
equation in quantum mechanics. Here, however, it is re- 
lated to the motion of Cooper pairs. Due to the diffusive 
nature of Cooperons, the mass entering the conventional 
Schrodinger equation is replaced by the inverse of the dif- 
fusion constant \/D. The solutions of this equation de- 
pend on the external magnetic fleld and on the geometry 
of the system, as the equation needs to be supplemented 
with appropriate boundary conditions. For the boundary 
to an insulator or vacuum, the following condition should 
be chosen 



n • (V - 2ieA)(j) - 0, 



(2) 



where n is the vector normal to the boundary. This 

condition corresponds to the requirement of zero super- 
current through the boundary. 

The information about the eigenfunctions (f)n and cor- 
responding eigenvalues q;„ for a certain geometry is suf- 
flcient in order to obtain the fluctuation corrections. 



We now present the results for the fluctuation correc- 
tions (in this section and the rest of this manuscript, we 
set h = 1). We write the total correction as the sum of 
three parts. 



5(T = dados + 5a an + 5a sc, 



(3) 



corresponding to the classification in the Usadel equation 
scheme. The density of states contribution takes the form 

5ados (4) 

= 2De2 f^J^Pnn [B'ReS'^ImLn - Blm {£'^1^)] . 

n 

This correction originates from the suppression of the 
quasiparticle density of states near the Fermi surface. 
The anomalous Maki-Thompson correction reads 



6a„ 



2De 



I 



TT-yZPnn— ImL^Imf^ . (5) 
2n ^ an 



This correction may be interpreted as a resonantly en- 
hanced interference effect in the Cooper channel. 

The correction induced by the fluctuating supercurrent 
is conveniently written as the sum of three terms 



(6) 



where 



J 277 f-f a™ - a„ 



and 



dlO v ^ dnrriB' 

2lT ^-^ am — ar. 

nra 



Im (f „ - f „ 



xlm (L„ - Lm) , 

xRe {Ln -Lm + - £m) LnL*J . 



(7) 



(8) 



(9) 



We introduced the retarded fluctuation propagator in 
equilibrium 



Ln{uj) = 



1 



f„(w)' 



(10) 



where 



T 

(w) = In ^ -h V' , 2 



1 2an+uo' 
2^ 4wT ^ ' 



and ip denotes the digamma function. B {oj) = 
coth (w/2T) is the bosonic equilibrium distribution func- 
tion. The prime in the above set of formulas denotes 
a derivative with respect to frequency, /'(w) = duf{oj). 
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We further introduced the following matrix elements in 

the basis of eigcnfunctions: 

Pnm (r) = 0; (r) (j)m (r) (12) 

dnm (r) = (C (r) [V - 2ieA(r)]</)„ (r)) (e • r„„) (13) 

where E is the unit vector in the direction of the exter- 
nal electric field. For the derivation, we employed the 
relation = —dmn- This relation holds true for the 
longitudinal conductivity which we study here, but not 
for the transverse component. ^'^ 

We note that the index n, used iii general to enumerate 
the eigensystem of Eq. (1), might in fact be a multi-index 
with several components. It is also possible that n does 
not enumerate a discrete set, but rather a continuum. 
In that case, the sum over n has to be replaced by the 
corresponding integral. 

An important remark is in order here. The anoma- 
lous Maki-Thompson correction diverges in the absence 
of a magnetic field, as then ao 0. The correction 
may be regularized by introducing a finite dephasing rate 
l/r^.^ Dephasing can be provided by magnetic impuri- 
ties, electron-electron or electron-phonon collisions. For 
low temperatures, electron-electron collisions dominate. 
Outside the region of strong fluctuations, one can con- 
sider the dephasing rate as energy-independent and equal 
to the sum of rates due to the Coulomb^^ and Cooper 
channels. ^^'^'^ In our study, we will treat l/r^ as a phe- 
nomenological parameter; it may be be introduced into 
the theory by replacing a„ — )• a„ -|- l/2r0 in the formulas 
for 6a and £n given above. 

In principle, the equations presented above can be used 
to obtain the fluctuation corrections to conductivity for 
any sample along an unconfined direction, by solving the 
eigenvalue problem (1) for the given geometry. In the 
next section, we briefly list a number of cases, for which 
these formulas can be applied. 

3. Eigenvalues a„ for different geometries 

In a bulk sample, in the absence of a magnetic field, the 
eigenvalue equation (1) can be solved by Fourier transfor- 
mation due to translational invariance, resulting in the 
following eigcnfunctions and eigenvalues: 

fAr)=e'^-, a, = \Dq\ (14) 

Films and nanowires have a reduced dimensionality. 
For the unconfined directions, it is useful to introduce 
a continuous Fourier transformation, while in the trans- 
verse direction(s). modes are quantized. This remains 
true if a parallel field is applied, because it is possible 
to introduce a vector potential that depends only on the 
transverse coordinate(s). The eigenvalues 

Ctqn = \Dq^ + a^n (15) 



can be written as a sum of a continuous and a discrete 
component. 

If the transverse direction is small in extent com- 
pared to the superconducting coherence length often 
only the lowest transverse mode is relevant. The low- 
est eigenvalue aj^o then plays the role of a pair-breaking 
parameter. ^'^^ For example, for a film in a parallel mag- 
netic field one finds 

a^o = \D{eBfd^, (16) 

where d is the thickness of the film. 

For a film with a perpendicular field, the situation is 
different from the one discussed above, since the vector 
potential A cannot be chosen to be translationally in- 
variant within the plane. The eigenvalues are degenerate 
Landau levels 

an = 2eDB + ^) - (17) 

which should be supplemented with suitably chosen 
eigcnfunctions.^^ For this case, the fluctuation conduc- 
tivity has been discussed with the help of the Usadel 
equation approach in Ref. 20. 

In this paper, we show that in a tilted field with both 
perpendicular and parallel field components, the relevant 
eigenvalues a„ can be written as a sum of the eigenvalues 
for the perpendicular and parallel magnetic field cases. 
We will discuss the fiuctuation corrections for this case 
in detail. 

Another interesting quasi one-dimensional system is 
a cylindric shell, i.e., a nanowire with annular cross- 
section.^'^ Here, a parallel magnetic field also gives rise 
to a flux threading the cylinder. The dependence of the 
lowest eigenvalue a±o on the flux is then (in the limit of 
vanishing thickness) periodic with the superconducting 
flux quantum (p^ = h/2e, as can be seen from solving the 
eigenproblem (1). This special case will be discussed in 
a separate publication.^^ 

B. Film in a tilted magnetic field 

We will now discuss in detail the fluctuation correc- 
tions for a thin film in a tilted magnetic field. We con- 
sider a thin amorphous superconducting film of thick- 
ness d penetrated by a magnetic field at an angle 6, 
< 6 < tt/2, measured between the field lines and the 
sample. We study the dirty limit, i.e., Tcqt <C 1, where 
Tco is the critical temperature of the superconductor and 
T is the elastic scattering time of the electrons. We choose 
coordinates so that the film lies in the x — y place, and 
the magnetic field B can be written as (see Fig. 1) 

B = Bsm9z + Bcosey. (18) 

Here, z and y are unit vectors in the z and y direction, 

respectively. We will sometimes use the notation By = 
Bcos{9) and Bj_ = Bsm{0). 
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Figure 1. A sketch of the system under study. A supercon- 
ducting film of thickness d is penetrated by a magnetic field 
of magnitude B at an angle 6. 



The film is assumed to be sufficiently thin so that the 
condition d £_o is fulfilled, where is the supercon- 
ducting coherence length at zero temperature. In this 
limit, the film can be considered as two-dimensional as far 
as its superconducting properties are concerned, whereas 
the electron motion is assumed to be three-dimensional. 
In the following we neglect the destructive effect on su- 
perconductivity caused by the direct coupling of the mag- 
netic field to the magnetic moment of the electrons. It is 
known that this is a good approximation for perpendicu- 
lar magnetic fields and weakly disordered films, epT 3> 1, 
where ep is the Fermi energy. ^ For parallel magnetic fields 
there is a minimum thickness ddog = ^o/^ft below which 
paramagnetic effects start to dominate. This is known as 
the Clogston limit, and we will assume that it is not 
reached, i.e., that d 3> ddog- 



1. Phase diagram, 

The experimental phase diagram of the film is three- 
dimensional. It is spanned by the temperature and the 
magnetic field, which is further characterized by its mag- 
nitude B and direction 9. A cut of the mean field phase 
diagram for zero temperature is displayed in Fig. 2. The 
phase boundary of the superconducting film is deter- 
mined by the pair-breaking parameter a, which quan- 
tifies the effectiveness of the magnetic field to suppress 
superconductivity.^ In the considered case of a tilted 
magnetic field, the pair-breaking parameter is the sum 
of the contributions due to the parallel and perpendicu- 
lar field components, a = an + a_\_, where 



a±^ — DeB sin ( 
Dd^ 



(19) 



6 



{eB cos ( 



Here, D = ^VpT is the electronic diffusion constant of 
the material, where vp is the Fermi velocity. 

The critical pair-breaking parameter ac, which sepa- 
rates the normal phase with a > from the supercon- 
ducting phase with a < ac, is temperature-dependent. 
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Figure 2. Phase diagram of the film in a tilted magnetic field 
at zero temperature. The diagram is parametrized by pa- 
rameters h and z, illustrated by lines of constant h (solid) 
and constant z (dashed). SC denotes the region of supercon- 
ductivity. Region I and II refer to the asymptotic regions in 
which we evaluate the corrections to conductivity. They are 
separated by the line h = 2/(1 -\- z). The sharp angular de- 
pendence of the negative correction to conductivity discussed 
below is a result of the pinching of region II near B^w ■ 



It is implicitly defined by the equation 

T 

. Tco 



In 



+ V - + 



0, (20) 



where tp (x) is the digamma function.^'* This equation 
has no solution for temperatures T > Tcq, for which the 
system is a normal metal. By applying the asympotic 
expansion 'ijj (x) « Ina; for large x, one finds that for zero 
temperature aco — ^'f^i where ln7 « 0.577 is the Euler 
constant. 

Once the solution ac (T) is known for arbitrary tem- 
peratures, the phase boundary can be constructed. It is 
a two-dimensional surface in the three-dimensional phase 
diagram spanned by the perpendicular and parallel com- 
ponents of the magnetic field B± and and by tempera- 
ture. The phase boundary is determined by the equation 



aj_ + an = ac (T) . 



(21) 



In order to find the critical field as a function of tem- 
perature at a fixed angle 6, for example, one should 
insert the expressions for a± and ay of Eq. (19) into 
Eq. (21) and obtains a quadratic equation for the critical 
field strength Bc{0,T). For a cut at constant tempera- 
ture T < Tco, it is convenient to present Eq. (19) in the 
form^-^'*'^ 



B, 



{9,T) cose 
Bc\\{T) 



Be {0, T) sin t 
Bc±{T) 



-1, 



(22) 



where Bc± and Be 
critical fields for 9 = 



are the (temperature-dependent) 
7r/2 and 9 = 0, respectively. They 
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can be obtained by setting q:||,_l = olc{T) and resolving 
for B. Eq. (22) describes a parabolic phase boundary in 
the (Bj^ — _Bj[)-plane. 

So far, the effect of a finite dephasing time was not 
included. As pointed out before, it may be accounted for 
by a shift in the eigenvalues a„ — >■ Q!„ + 1/2t0. Therefore, 
the condition for the mean-field transition can be written 
as a = ac{T), where ac{T) = adT) — 1/2t^ is modified 
due to the presence of dephasing effects. If the dephasing 
time is weakly magnetic-field dependent, then its main 
effect on Eq. (22) is to renormalize the critical fields B^^^ 
and Bc±. 

We note that Bc± [T = 0) coincides with the nucleation 
critical field Bc2 ■ For low temperatures, T — >■ 0, the two 
critical fields B^u and i?c_L are related by 



' -J =4.16x5^, (23) 

I V ttco d 



Se±(T = 0) 

where is the zero-temperature coherence length.^ 



2. Parameterization for the vicinity of the quantum critical 

line 

When formulating the results for the film in a tilted 

magnetic field below, we will specifically discuss the vicin- 
ity of the quantum critical line in the phase diagram, i.e., 
we concentrate on low temperatures, t = T/Tco <^ 1. For 
the fluctuation conductivity, the regime of small tem- 
peratures is particularly interesting. This regime dis- 
plays the phenomenon of the non-monotonic magnetore- 
sistance (NM), as was first shown for the perpendicular 
magnetic field in Ref. 11 and for the parallel magnetic 
field in Ref. 12. Here, we will discuss this phenomenon 
for magnetic fields tilted at arbitrary angles 6. 

For a fixed temperature, the phase boundary has the 
shape of a parabola, cf. Eq. (22). In order to present 
the results for the fiuctuation corrections, we introduce 
a paxametrization of the region close to this critical line, 
which will be described in the following. 

First, choose an arbitrary angle 9 between and 90 
degrees. Eq. (20) and (22) determine the critical field 
strength Be {9), at which, for a given angle, the system 
undergoes the phase transition. The strength of the mag- 
netic field can then be measured by the relative distance 
h to the phase boundary: 

B-B, (9) 



h = 



(24) 



Be (9) ' 

with h = corresponding to a point on the critical line 
and h > corresponding to a point in the normal region 
of the phase diagram. 

In order to parametrize the angle 9, it turns out to 
be useful to introduce another dimensionless number z, 
which is defined as the ratio between the two projected 
pair-breaking parameters 

ceil eBdP cos^ 9 



As can be seen from the definition, z is directly related 
to the angle, with z = corresponding to perpendicular 
field and z = oo corresponding to parallel field. The 
parameters h and z can be used instead of B and 9 in 
order to define a point in the phase diagram. In fact, 
one can consider h and z as a new curvilinear coordinate 
system of the phase space that is aligned along the phase 
boundary of the system, as depicted in Fig. 2. 

Whenever the vicinity of the transition line in the 
[B±_ — i?||)-plane for a fixed temperature T <C T^q is con- 
sidered, the B dependence of z may be neglected and one 
may approximate z w z{B = Be {9,T) ,9). This quan- 
tity can be determined experimentally (without explicit 
reference to the thickness d) : 



.(5 = s.(.,T),.) = ^^J;^-i. 



(26) 



ziB,9) = 



aj_ 



sin 9 



(25) 



For almost parallel magnetic field, ^ <C 1, and for T ^ 0, 
one finds 

Be±(T = 0)1 dl 

z=—^ ^-= 0.24 x . (27) 

B,||(T = 0) 9 " Cod ^ ^ 



3. Results for the fluctuation conductivity 

The kinetic equation approach employed in this paper 
leads to a rather natural classification of the three dis- 
tinct contributions to fluctuation conductivity according 
to the underlying physical mechanisms. We distinguish 
the density of states contribution, Sudos, the anoma- 
lous Maki-Thompson term doani and the contribution 
of the fluctuating supercurrent 5(Jse, as discussed in the 
introduction. It is worth mentioning that this classiflca- 
tion differs from the conventional diagrammatic scheme. 
For the parallel magnetic fleld case, the precise corre- 
spondence between the two formalisms is worked out in 
Sec. IV, and illustrated in Tab. I. 

We will now state results for the fluctuation conduc- 
tivity of the film in the vicinity of the quantum critical 
line, i.e., for T <C Teo (f <C 1). Let us stress again that 
the general results displayed in Sec. II A cover the entire 
normal part of the phase diagram [with the exception of 
the region of strong fluctuations very close to the tran- 
sition]. Here, we focus on the low-temperature regime 
since it displays the interesting phenomenon of the NM. 
Whenever possible, we will discuss the origin of the differ- 
ent corrections according to the classiflcation into density 
of states, anomalous Maki-Thompson and supercurrent 
contributions. 

When formulating our results we make use of the pa- 
rameterization of the phase diagram in terms of the pa- 
rameters h, z and T introduced above, see Fig. 2. For 
the low temperature regime, it can be expected that the 
presence of in the fluctuation propagator mainly leads 
to a shift in the critical line in the (Bj^ — B\\) plane. We 
assume that this shift has already been performed. At 
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the same time, we neglect in the Cooperon, because at 
low T and in the vicinity of the critical line, the Cooperon 
is not singular. This is why the presented formulas will 
not contain any explicit reference to t^. 

Comparatively simple semi-analytical expressions can 
be found in two regimes. In region I, defined by the rela- 
tion h -C , Landau level quantization of the Cooper 
pair motion is crucial. Indeed, the dominant contribu- 
tion to fluctuation conductivity in this regime originates 
from fluctuations of the lowest Landau level, since these 
fluctuations become singular at the transition. In regime 
II, for which h ^ the spectrum may be approxi- 

mated by a continuum for the purpose of the calculation. 
The reason is that in this regime either, for small z, the 
distance to the transition line is comparatively large and 
the fluctuations of all levels are non-singular or, for large 
z, the magnetic field is almost parallel and the distance 
between adjacent Landau levels becomes very small. The 
two regions are displayed in the diagram of Fig. 2. 

With the only exception of very small angles 6, the 
system is in region II of the phase diagram when ap- 
proaching the transition line, compare Fig. 2. In turn, 
for large magnetic fields far from the transition li ^ 2, 
the system is in region II independent of the value of 
z. The cross-over angle between the two regimes near 
criticality can be estimated from Eq. (27): Region II is 
reached only for very small angles 



<C 0.12 X '^h. 



(28) 



As an illustration, for a film of thickness d = 0.3 ^o, fairly 
close to the transition h = 0.1, the cross-over occurs at 

an angle 6 of about 0.21°. 

In the following, we will discuss the two regions sep- 
arately, starting from region I. After stating the results, 
we will provide a qualitative discussion of the behavior 
in the two regions. 

Region I: As a special example, region I contains the 
case of a strictly perpendicular magnetic field, z = 0, 
in the the vicinity of the transition. This case has first 
been treated in Ref. 11 (see also Ref. 20). Our results 
for region I can be viewed as a generalization of these 
previous results to non-perpendicular angles. 

The general formulas stated in Sec. II A involve an 
integration over an internal frequency and the summa- 
tion over Landau level indices. As was already noticed 
in Ref. 11, the most singular contribution in the vicin- 
ity of the transition stems from the lowest Landau lev- 
els (LL) only. Correspondingly, we consider the contri- 
butions due to the singular LL and due to the higher 
Landau levels (HL) separately. It turns out that for 
the HL a continuum approximation is suflicient. Fur- 
thermore, for both LL and HL contributions we perform 
a separation into a thermal correction (T), which van- 
ishes for T ^ 0, and a quantum correction (0), which 
is temperature-independent and thereby persists even in 
the limit T — )■ 0. As a result, the corrections to conduc- 



tivity may be presented in the following form: 
6ai = Sao,LL + SaT,LL + Sao,HL- 



(29) 



Here, the thermal contribution Sar^LL reads as follows 

2 

SaT,LL = ^ (r) + PI^ (r)) , (30) 



where 



/„=lnr---^(r), (31) 
2r 

If> = r^' (0 - ^ - 1, (32) 



and we have abbreviated r = tAt. This contribution is 
very similar in structure to the result for the perpen- 
dicular magnetic fleld case derived in Ref. 11. It differs 
mainly in two respects. First, we omitted the term In/i 
from la as it does not vanish as T — >^ and is therefore 
part of the quantum contribution Sao,LL to be discussed 
below, cf. Eq. (36). Second, the prefactors a and (3 are 
now z-dependent, i.e., they depend on the angle 0, 



1 



5 + 2z 



a 



3 + z 1 + z' 



+ 



1 



3 + z 1 + z 



(33) 



In the case of a perpendicular magnetic fleld, the coef- 
ficients reduce to the previously derived a = —2/3 and 
/? = 8/3.^^'^° Interestingly, however, in approaching the 

transition, /i — > 0, we find that ^ctt.ll ~ '^^^S-j^, mean- 
ing that the z-dependence of the general formula drops 
out in this limit. When approaching the transition at 
any finite temperature, S(Tt.ll eventually becomes the 
dominant contribution. It then resembles the well-known 
Aslamazov-Larkin fiuctuation correction. ^^'^^ 
The quantum contribution due to the LL reads 



0,LL 



1+z 1+h V z+1 



(34) 



where Li is the logarithmic integral function.^* The first 

term stems from the density of states correction Sados-, 
and the second part is due to the supercurrent correction 
8asc- The function G {h, a) is defined by the integral 

G{h,a) = \ r dx , , f-^ - — ^ (35) 

^ ' 2 Jxj^h x{x + a)\\xi.x \a.x + a)^ ' 

0, the contribution Jctq ll can be seen to 



For h 
reduce to 



LL 



s^aln 



h 



(/i^O), 



(36) 



which corresponds to the quantum term in the formula 

found by Galitski and Larkin.^'^ 

The last contribution in Eq. (29), Sao,HL, was omit- 
ted in Ref. 11, and is obtained by considering the higher 
Landau levels. For this term, one may use the contin- 
uum approximation for the sum over Landau levels. This 
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formally corresponds to the limit z — > 0, i.e.. this term 
is only very weakly z-dependent. In addition, it is not 
singular when approaching the transition. Formally, the 
sum over higher Landau levels is very weakly (doubly 
logarithmically) divergent, so that it becomes necessary 
to introduce a high-energy cut-ofF A and to take into ac- 
count only modes with a„ < A. As our theory is based 
on the diffusion approximation, the cut-off can be cho- 
sen to be of the order of the transport scattering rate, 
A~ 1/r. 

The quantum correction from higher Landau levels can 
then be written as the sum of two integrals: 



c- (dos) 



(37) 



27r2 



dy 



K 



dx 



{a + x + y) \B.{a + x + y) 



6a, 



(scl) 



e 

27r2 



0,HL 
1 



(38) 



dy 



K 



dx 



(a -I- a; + y^ \t? {a^ x ^ y) 

where we abbreviated a = (1 + /i)^, and K is a. dimen- 
sionlcss cutoff Given by K = t^^. These integrals can 
be explicitely solved in terms of the logarithmic integral 
Li, see Eq. (85) in Sec. V. 

As a final remark concerning the correction in region I, 
the thermal contribution originating from higher Landau 
levels, SaT,HL, has been omitted from formula (29) as it 
is regular and small. 

Region II: We now turn to the discussion of region II, 
where Landau levels are so close that they can entirely 
be treated in the continuum approximation. In that case, 
the special significance of the lowest Landau level is lost. 
As nicntionecl earlier, the limit of taking a continuous 
spectrum corresponds to the limit 9-^0. The results 
for a film in a strictly parallel field have been found pre- 
viously in the diagrammatic technique, and we agree 
with these results. For a more detailed discussion of the 
comparison, we refer to Sec. IV. 

We separate the total correction to conductivity in re- 
gion II into thermal and quantum contributions: 



Sar 



6aT, 



HL 



5ao, 



HL- 



(39) 



The expression for the quantum contribution Sao^HL is 
the same in both regimes I and II and has already been 
stated above in Eq. (39). The thermal part is domi- 
nated only by the supercurrent correction Sai^''^ in the 
classification introduced in Sec. II A [in the conventional 
classification, it originates from the Aslamazov-Larkin di- 
agram]. It can be written in the form^^ 



Sax. 



HL 



F (ry) . 1] 



ac (T) 



T 



(40) 



For sufficiently small h, we can approximate ry = 27rr. 
The function F in the previous formula is defined by 



Fin) 



ij/ (x) dx 

X 



1 

47r 



1 
16?7 



(41) 



It may be evaluated numerically, the asymptotic ex- 
pansion for large and small values of the argument gives 



for ry ^ 1 
for ?y <C 1. 



(42) 



We now turn to a more qualitative discussion. 



4- Qualitative discussion: Film in a tilted magnetic field 

Region I: The two dominant corrections in region I 
are SaT,LL and Sgq ^l. They are of different sign when 
approaching the transition. For small /i <C t, in the so- 
called thermal regime, SaT,LL dominates. All the cor- 
rections Sados, Saan and ^^sc contribute to SaT,LL and 
the net result is a positive correction to conductivity as 
naively expected when approaching the superconducting 
state. However, the presence of the term 5ao,LL, which is 
negative and dominates in the so-called quantum regime 
t, leads to a local minimum in 5a as a function of h 
at the cross-over between the thermal and the quantum 
regime. The result is the NM.^^'^^ The physical origin is 
the negative density of states correction, which competes 
with the positive supercurrent correction in the quantum 
regime h ^ t, but is numerically larger. The anoma- 
lous Maki-Thompson term is ineffective in this regime. 
The physical mechanism underlying the phenomenon of 
the NM in region I can therefore be expressed in simple 
terms. Cooper pairs form, but are comparatively immo- 
bile as a consequence of their quantized spectrum (Lan- 
dau level quantization). The decrease in conductivity 
due to a reduction of the density of states of quasiparti- 
cles may therefore overcome the increase in conductivity 
caused by the fluctuation supercurrent carried by these 
Cooper pairs. 

Crossover between regions I and II: The importance 
of the term Sao^HL fies in the fact that, irrespective of 
the precise choice of the cut-off, it is negative and weakly 
angular dependent. Deep in region I, the phenomenon of 
the NM is largely determined by the interplay of 6ao,LL 
and SaT,LL as discussed above, and Sao,HL is of minor 
importance. As one approaches the cross-over regime 
between regions I and II at small angles, however, the 
negative correction due to Sao^LL diminishes, as it is pro- 
portional to a, and Sa^^uL becomes more relevant. As 
will be discussed in more detail below, in region II the 
NM still exists only due to the existence of the negative 
contribution Sao,HL- For the accurate description of the 
cross-over regime itself, a numerical evaluation is neces- 
sary and results are displayed in Fig. 3. Let us remind 
at this point that the negative correction 5ao,HL stems 
from the density of states and supercurrent corrections 
5a dos and 5a sc- Again, the negative density of states 
contribution dominates. 

Another interesting observation can be made in the 
thermal regime /i <C i, which can be reached in both 
regions by approaching the transition. It is characteristic 
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for this regime that the divergent thermal contributions, 
Sct^ll and SaT,HL, dominate. As mentioned before, very 
close to the transition, the system is in region I for almost 
all angles. The cross-over to region II only happens very 
close to 6* = 0. It is thus interesting to note that for h ^ 
t, the asymptotic expansion of Sut.ll contributing to 
5(ji is while the asymptotic expansion of SaT,HL, 

which contributes to Sau, is Thus, while crossing 

over from region I to region II in the thermal regime, the 
fluctuation correction drops to about 1/8 of its value for 
small h. 

Region II: As mentioned above, Sao.HL is regular at 
the transition, and slowly decreases when moving to- 
wards the normal regime. On the other hand, the ther- 
mal contribution 6aT,HL diverges when approaching the 
transition. As demonstrated in Ref. 12, the interplay of 
these two contributions also results in a NM. From the 
preceding discussion it is clear that the non-monotonic 
behavior of the magnetoresistance has a different origin 
for parallel and perpendicular magnetic fields. For paral- 
lel fields, the negative correction comes from Sao HL- For 
perpendicular fields, it originates from Sao,LL- In both 
cases, however, these negative corrections stem from a 
competition of density of states and supercurrent terms, 
for which the negative Sados dominates over the positive 

S(Tsc- 

This concludes the discussion of the main results. In 
the following sections, details of the derivation will be 
presented. 




0.00 0.01 0.02 0.03 

h = 

Be 

Figure 3. The correction to conductivity of a disordered film 
at very low temperatures (T = 0.005 T^o, d/^o = 0.3). The 
transition is approached by tuning the magnitude of the mag- 
netic field at a fixed angle 9. When one tilts the angle starting 
from the perpendicular magnetic field case, the dominant ef- 
fect is a reduction of the negative correction to conductivity. 
When tilting further towards parallel fields a second effect 
sets in, namely a sudden decrease of the thermal contribu- 
tion. This effect becomes visible mainly for small h. 



III. GENERAL FORMALISM 

In this section, we describe the formalism underlying 
the results presented in Sec. II A. We derive a quasiclas- 
sical kinetic equation, the so-called Usadel equation. The 
form of this equation is slightly nonstandard in order to 
allow for the inclusion of fluctuations outside the super- 
conducting regime. The formalism used here was intro- 
duced in Ref. 20 for a film in a perpendicular magnetic 
field. We generalize it here in order to be able to treat 
parallel field components as well. The main difference to 
Ref. 20 is that a more general set of eigenfunctions and 
eigenvalues is considered. Otherwise, the derivation is 
analogous to the one presented in Ref. 20. We describe 
the main steps here in order to coin the notation and to 
make the paper self-contained. Previous related works 
include Ref. 27, in which the Usadel equation was used 
for the calculation of fluctuation corrections close to Tco 
and Refs. 44 and 21, where the Keldysh nonlinear sigma- 
model was employed for the calculation and again, only 
the vicinity of Tco was studied. The latter two works are 
based on the Keldysh sigma-model approach for super- 
conductors presented in Ref. 45; the connection between 
the different formalisms is that the Usadel equation is the 
saddle point equation of the sigma-model. 



A. Microscopic model 

We start from the Keldysh action for electrons with 
BCS short-range interaction. The interaction in the 
Cooper channel is already decoupled by means of a 
Hubbard-Stratonovich transformation, leading to the ac- 
tion 




+A*xat + Ax;x: + ^A*A. (43) 

Here, xt(^) ^^^d Xii^) are Grassmann flelds describing 
the spin-up and spin-down components of the electrons 
at the space-time point x — {r,t), h is the single-particle 
Hamiltonian, A(x) is the superconducting order parame- 
ter fleld, 1/ is the density of states of electrons per spin di- 
rection at the Fermi surface, and A > is a dimensionless 
parameter, which determines the strength of electron- 
electron attraction in the Cooper channel. The integral 
in time is along the closed Keldysh contour C*^ The sin- 
gle particle Hamiltonian 

h = —(V-ieA(r)f + U(r)+eLp(r) (44) 

2m 
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includes coupling to a vector potential A and an external 
scalar potential (p, as well as an impurity potential U. 
In this work, wc assume a constant potential gradient 
ip (r) = — E ■ r. where E is the electric field. 

We separate fields on the forward and backward branch 
of the Keldysh contour. We further introduce the Nambu 
(N) particle-hole space, and perform the Keldysh rota- 
tion. This means that we describe our system by the 
four-component fields * = (f/'i, ^'2)^' = (Xti-Xji)!^ 

and = (V^, '^'''|)a'. V4 = (Xt,:- -X-i-^)^, where K labels 
the Keldysh space. This transformation turns the action 
into 



-^tr(AtaiA). 



(45) 



Here and in the following, fj denote Pauli matrices in 
particle-hole space, while CTj denote Pauli matrices in 
Keldysh space. The trace is taken in the four-dimensional 
product space. Matrices in this product space are de- 
noted as A, H, etc.. The time-integration covers the real 
axis. The single-particle Hamiltontian H is obtain from 
the Hamiltonian h by replacing A by Afs . The 4x4 ma- 
trix A contains the classical (c) and quantum (q) compo- 
nents of the order parameter field A as A = Aqcto-I-Aicti, 
where Aj = Ajf+ — A|f_, t± = (f^ ± «fy)/2 and 

A=(A„A,)^.20.46 

In order to perform the quasiclassical approximation, 
it is useful to define the electronic Green's function 



iGa (x, x ) = — — — — 



,iS[*,A] 



(46) 



where the average is performed over the electronic 
degrees of freedom only, keeping the field A fixed. 
In order to recover the full Green's function G one 
needs to average as G = J [dA] Gac'^^^'I^J with 
respect to the Ginzburg-Landau action Sql [A] = 

— i In J [(i^"] . In Keldysh space, the matrix 

Green's function has the typical triangular form 



G 



K 



G^ G 
G^ 



(47) 



It is important to note that in general Ga does not have 
the same structure. If Ag ^ 0, the element in the lower 
left corner, sometimes referred to as G^, is not equal to 
zero. For the linear response calculation of the fluctua- 
tion conductivity in the Gaussian approximation, how- 
ever, it turns out to be sufficient to put G^ = as it 
is proportional to higher powers of the order parameter 
field. This has been shown in Ref. 20. We will therefore 
work with the triangular Keldysh structure for Ga. 

In the presence of impurities, the Green's function 
G needs to be averaged over an ensemble of disorder 
configurations. Physical quantities can be calculated 



with the help of the disorder averaged Green's func- 
tion (Ga)^^^^. For films with a dimensionless conduc- 
tance g ^ 1, which we consider, it is legitimate to av- 

eragc as (G),^^ = / [dK] (Ga),,,^''^^"^^^-, i.e., sep- 
arately for Ga and Sgl- The quasiclassical approxima- 
tion may now be introduced for the disorder averaged 
Green's function {G^^ns at a fixed order parameter con- 
figuration. 



B. Quasiclassical approximation and the Usadel 
equation 

We assume that ep is the largest energy scale in the 
problem, meaning that superconductivity, scattering and 
external fields only affect the system close to the Fermi 



surface: A, t 



<^ ep. Thus, one may use the so- 



called quasiclassical approximation.^''' To this end one 
transforms to Wigner coordinates and defines the quasi- 
classical Green's function by integrating over the distance 
to the Fermi surface 



5n (r,i,t') 



.,{r,p,t,t'). 



(48) 



Here, n is a unit vector pointing in the direction of p, 
2 

and = 2m ~ quasiclassical Green's function g„ 

obeys the Eilenberger equation. 

In the limit of dirty superconducitivity, Tcqt <C 1, fre- 
quent scattering washes out the angular dependence of 
the Green's function, which allows to work with the sim- 
pler function g = J dxig^- The appropriate equation in 
this limit is the Usadel equation^^'*^ 

DdA (gdAg) - {hdu g}+i[K-^,g]= 0- (49) 

Here, D = VpT/3 is the electronic diffusion coefficient 
for diffusion in three dimensions and ip is the external 
potential, which is considered as a matrix in Nambu and 
Keldysh space. We defined the derivative 9a5 = — 
ie [Afs, g]. Also note that one should read the appearing 
anti-commutator in an operator sense, {fsdt, g} (^1,^2) = 
T39tig{ti,t2) — dt2g{ti,t2)T3. In this equation and the 
following, we consider g, A and (p as matrices in Keldysh 
and Nambu space and in time. Thus a mTiltiplication 
implies the matrix product of the 4x4 matrix as well 
as a convolution in time. The Usadel equation has to be 
supplemented with the constraint g^ = 1. The current in 
the sample can be expressed through g as^^ 



(t3(Ti (^gVg 



(50) 



This procedure is an extension of the usual quasiclassi- 
cal formalism in the Keldysh technique; here, the formal- 
ism is applied without employing the mean-field approx- 
imation for the field A. This allows us to treat the elec- 
tronic system in the quasiclassical approximation while 
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fluctuations of the order parameter field A can still be 
taken into account. 

In the rest of this section, we consider the normal side 
of the transition, where fluctuations of the order param- 
eter fleld are taken into account in the Gaussian approxi- 
mation. We therefore linearize the Usadel equation using 
the solution A = as a starting point. This allows to 
obtain an expression for the current as a functional of 
the field A, which can then be average using the corre- 
lation function for A. It is possible to obtain the corre- 
sponding fiuctuation propagator from the quasiclassical 
Green's function itself. 

In analyzing fiuctuations around the normal state of 
the metal, we parametrize the Green's function as in Eq. 
(47), using 



1 - ff*/2 

r 



-i + r//2 



-1 + fr/2 -J 
-I* 1 - r//2 



(51) 
(52) 



while for the Keldysh component we set g = g h — hg . 
In these equations g^ , , g^, h, and / are all functions 
of one spatial and two time coordinates, = g^{r, t, t') 
etc., and the product implies a convolution in time. The 
paramctrization is consistent with the non-linear con- 
straint = i up to third order in /. The functions /,/*, 
/, /* are considered to be independent from each other; 
in fact, it turns out that /* is the conjugate of / only if 
Aq = 0. Inserting (51) and (52) into the Usadel equation 
(49), we obtain a linear equation for f,C~^f = V, where 



Vt„t2 (r) = 2i[Ae(r, ti)St,-t, + he{r, h - ^2) Ag(r, is)], 

(53) 

the operator is defined as 

C-^ =Ddi-{dt,-dt,), dA = V-2ieA{r), (54) 

and he will be introduced in Eq. (55) below. For /*, /, /* 
similar equations are obtained, which show that fti,t2 = 
-/t2,ti> ftut2 = ^liile in general /* ^ (/)*. ' 

From the kinetic equation for h, to lowest order in A, 

one finds that the only deviation from the equilibrium 
distribution function is a shift due to the local potentiaP° 

h (r, e) = diag (ft.e(r, e), hh{r, e)) , (55) 
K/hir,e)^n{eTeip{r)). (56) 

Here, 'H(e) = tanh(e/2r) is the fermionic equilibrium 
distribution function. 

In order to solve the linearized Usadel equation for /, 
we introduce the Fourier transform 



ftiM — I 



ieiti+ie2t2 



(57) 



where J^ = J de/2'K. The same convention is used for the 

independent variables /*, /, /*, and for V. In a second 
step, we expand / = En /nV'n, /* = En /nV'n^ etc., and 



use that the spatial part of the operator C is diagonal 
in the basis of eigenfunctions i/'n of Eq. (1). It gives 

/n;e„e. = C„;e,+., j dv rni^) V;„,,(r), (58) 

where we introduced the so-called Cooperon C„ 

C-^, =ie- 2a„ - (59) 

including a phenomenological dephasing rate as discussed 
in Sec. II A. 

We now turn to the fiuctuation propagator £, which is 
directly related to the correlation function of the field A 



{Ai{xi)A*{x2)) = ^^ij{xi,X2). 



(60) 



The indices i, j label the components of the vector A 
defined in Eq. (45). Conveniently, the fiuctuation prop- 
agator can be found from the quasiclassical Green's 
function. ^° We restrict ourselves to the Gaussian approx- 
imation for the fiuctuations of the order parameter field, 
i.e., we find the correlator by approximating Sol by its 
expansion to second order in A. The resulting term 
contains the Greens function Ga evaluated at coincid- 
ing space points, and can thus be expressed through the 
quasiclassical Green's function g. One finds 

(£-1) (r, r', t, t') = ^"jX^^^^^i^-a (r, t, t) 



ai] 5ir-r')6{t-t'), (61) 



2u 



For the first term on the right hand side the identification 

(7c = 1 and (Tq = (7i is used. The fiuctuation propagator 
can thus be found from the solution of the Usadel equa- 
tion. It is convenient to work in the basis of eigenfunc- 
tions of (1) and to present the fiuctuation propagator 
in the form £(r,r') = Y.„„tpn{'r)Cnm'^Pmi^')- To first 
order in the electric field, one finds the following matrix 
elements 



2ujC^^ (w) = -I (A,,„ (uj) A;„ i^)) 
2vC,t^ (w) = -I (A,,„ {uj) A;„ {^)) 



(62) 



and 



= 13 {Ln — L*^) Snm 

-2BEr„„(£:;i„i„-0:im)- (64) 

Furthermore, Tnm are the matrix elements of the position 
operator 



j drril)l (r) 



"^m (r) . 



(65) 
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and B (oj) = coth {u /2T) is the bosonic equilibrium dis- 
tribution function. L„ is defined in Eqs. (10) and (11). 
In Eq. (11), an should be replaced by a„ + 1/2t^ in or- 
der to account for the finite dephasing rate introduced in 
Eq. (59). 

In the next step, we insert the parametrization (51), 
(52) into formula (50). From the normal-metal solution, 
for which = —g^ = fa, one immediately finds the 
Drudc expression for the current j*^"^ = 2z/e^Z?E. Taking 
into account fluctuations up to second order in /, one 
obtains several terms with distinct structure. They can 
be grouped as follows^° 

.(do,) ^ £^ + ^ (6g) 

= -e\uDE [ ^'(e) (/«'/,%). (67) 

J ee' 

-/«'V+/.'.-V+ (68) 

where we abbreviated V± = V =F 2ieA. In the result 
(68), we dropped terms that contain 'H'{e). These terms 
vanish after averaging over fluctuations. The naming of 
these three components refers to density of states (dos), 
anomalous Maki-Thompson (an) and super-current (sc), 
respectively. 

Eqs. (66), (67) and (68) still need to be averaged over 
the fluctuating order parameter fleld. To this end, one 
needs to express /, /*, etc. in terms of A and A* and 
the Cooperon, as exemplified Eq. (58). The averaging is 
then performed with the help of Eq. (60) . The integral 
over the frequency e can be performed analytically. This 
step is straight-forward in principle, but tedious, this is 
why we display some intermediate steps in App. B. The 
final results for the fiuctuation corrections are displayed 
in Eqs. (3) to (9) in Sec. II A. 

IV. COMPARISON TO THE DIAGRAMMATIC 
TECHNIQUE 

In this section, we summarize the most important as- 
pects of a comparison between the results for the fiuc- 
tuation conductivity obtained in this manuscript, see 
Sec. II A, and those obtained with the help of the tra- 
ditional diagrammatic technique in Ref. 12. The com- 
parison is performed for a certain class of pair-breaking 
transitions considered in Ref. 12, for which our results 
are also applicable. As discussed in Sec. II A, one- and 
two-dimensional systems in a parallel magnetic field can 
be treated on a similar footing. The relevant eigenfunc- 
tions of Eq. (1) are plane waves (/'q(r) = e''^'', where q is a 
one- or two-dimensional wave vector. The corresponding 
eigenvalues have the form aq = \Dq^ +a, where a is the 
so-called pair-breaking parameter. In App. A it is shown 
that the two sets of formulas for the fiuctuation conduc- 
tivity, one obtained from the diagrammatic calculation 



and the other one with the help of the Usadel equation 
approach, can be transformed into each other by means 
of algebraic manipulations [the ultraviolet cut-off in 2d 
and 3d requires a special discussion, see App. A]. This is 
not only relevant from the point of view of consistency, 
but it enables us to find the relation between the classifi- 
cation used in Ref. 20 as well as in this manuscript, and 
the traditional classification based on diagrams. 

In the diagrammatic technique, one deals with five core 
diagrams, displayed in the leftmost column of Tab. I. 
Most of these have topologically similar partner dia- 
grams, so that in total a standard set of 11 diagrams 
needs to be considered, see Ref. 1. In these diagrams, 
full lines represent quasiparticle propagators, while wavy 
lines represent fiuctuation propagators. Disorder is taken 
into account with the help of the so-called impurity cross- 
technique,^^ shaded boxes and triangles, for example, 
symbolize Cooperons. 

The diagrams with one horizontal fiuctuation prop- 
agator are referred to as density-of-states (DOS) dia- 
grams, while those with vertical fiuctuation propagators 
are called Maki-Thompson (MT) diagrams. Both of these 
diagrams may contain either two (2c) or three cooper- 
ons (3c). The final diagram is the so-called Aslamazov- 
Larkin (AL) diagram. It contains two fiuctuation propa- 
gators. We would like to stress again that this classifica- 
tion is different from the one used in this manuscript. The 
authors of Ref. 12 further distinguish terms for which the 
frequency integral contains the factor coth(u;/2T) from 
those that contain 1/ sinh^((x)/2T). These are denoted as 
ch and sh, respectively, in our notation they correspond 
to terms containing the bosonic distribution function B 
or its derivative B' . The terms ^c^osc/s' ^'^mtcH 

^'^Dos ch' ^'^MTch ^^^J differ in their overall coefficient 
and in Ref. 12 they are called A and B, respectively. This 
naming scheme is also used here to refer to the terms ob- 
tained in Ref. 12, and they are displayed in the second 
column of Tab. I next to the corresponding diagrams. We 
now briefiy describe the manipulations required to con- 
vert the equations stated in Ref. 12 to our results. More 
details can be found in Appendix A. The first step is to 
perform the integration over fermionic frequencies e for 
all terms originating from the DOS- and MT-diagrams. 
The result of this integration can be expressed in terms of 
the propagator L. After this step, the algebraic structure 
is similar to the corresponding terms in the technique 
used in this manuscript and first identifications can be 
made. Finally, it is necessary to integrate by parts in 
the radial momentum variable q. These manipulations 
lead from column two to column three in Tab. I, and re- 
arranging the resulting terms to obtain the results from 
the Usadel technique leads from column three to column 
four. 

In summary, we were able to show that for a certain 
class of pair-breaking transitions there is a one-to-one 
correspondence between the results obtained from the di- 
agrammatic perturbation theory and those obtained from 
the Usadel equation approach. The results of this com- 
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parison are illustrated in Tab. I. 



We thus find the eigenfunctions 



V. FILM IN A TILTED MAGNETIC FIELD 

As an application of the formalism developed above, 

wc calculate the fluctuation conductivity for a thin super- 
conducting film penetrated by a magnetic field at an ar- 
bitrary angle. We focus on the low-temperature regime, 
where one can follow in detail the angular dependence 
of the resistance peak. The main results are summarized 
and discussed in Sec. II B. Here, we present details of the 
calculation. 

First, one needs to solve Eq. (1) in order to find the 
eigenmodes and eigenvalues specific for the sample ge- 
ometry and the orientation of the magnetic field. With 
this knowledge, the matrix elements Pnm and d„m of 
Eqs. (12) and (13) can be found. When inserted into the 
general expressions, Eqs. (3) to (9), one can obtain a set 
of formulas that would in principle allow to evaluate the 
fiuctuation corrections everywhere in the metallic phase. 
We will specialize on temperatures T ^ Tco, however, 
where one finds comparatively simple results. In order 
to structure the calculation, it turns out to be useful to 
distinguish quantum contributions, that are independent 
of temperature, and persist for T — > 0, and thermal ones, 
which capture the temperature dependence, but vanish 
in the limit T — >■ 0. Two asymptotic regions in the phase 
diagrams are found where the evaluation of the sum over 
Landau levels can approximately be performed. As the 
lowest Landau level is of special importance in certain 
limits, we evaluate the contributions of the lowest Lan- 
dau level and of the higher Landau levels separately. 



A. Eigenvalue equation and matrix elements 



Xr. 



X — 



2eB± 



(72) 



where are the eigenfunctions of a one-dimensional 
quantum harmonic oscillator with frequency 2eDB± and 
mass D~^, and corresponding eigenvalues 



a„ = 2eDBj_ 



1 



^p2n2 2 



(73) 



The Landau levels are degenerate, as the eigenvalues do 
not depend on p. The perpendicular component of the 
magnetic field causes Landau level quantization, and the 
parallel field component adds a constant to the pair- 
breaking parameter that is quadratic in It will prove 
useful to present the eigenvalues in the form 



a„ = a± {2n + 1 + z) , 



(74) 



where z is introduced in Eq. (25). 

The matrix elements required for the evaluation of the 
general formulas Eqs. (3) to (9) can be found by using 
well-known relations for the eigenfunctions of a quantum 
harmonic oscillator. As the eigenvalues do not depend 
on p, the integration over this variable can immediately 
be performed. We find 



Pnp^mq 



EB± 
2-K 



{mSn 



-1-1, m 



nS. 



n,m-t-l 



)• 



(75) 
(76) 



For the second matrix element, we only calculated the 
vector component in the x-direction, the direction of the 
electric field. This is sufficient for the calculation of the 
longitudinal conductivity. 



Here, we solve the eigenvalue equation (1) for a thin 
film in a tilted magnetic film. The tilting angle is de- 
noted as 6, for an illustration see Fig. 1 . We describe the 
magnetic field by a vector potential in Landau gauge 

A={-B^y + B\\z,0,0), (69) 

where we abbreviated B± = BsinO, B\\ = BcosO. The 
eigenvalue equation (1) then reads 

-^{{9x + "^ieB^y - 2ieB^^zf + 9^ + dl) = acj). 

(70) 

Under the assumption that the layer is thin, only the 
lowest transverse mode in the z-direction is of interest, 
for which is almost constant as a function of z. We 
drop d1 and integrate over z from —8/2 to s/2: 

((5. - 2ieB^yf + dl) <^^{a- ^e^B^s^^ cj>. 

(71) 



B. Low temperature approximation, two 
asymptotic regimes 

In the low-temperature regime t <^ 1, the calcula- 
tions can be considerably simplified by working with the 
asymptotic form of the fluctuation propagator. The cal- 
culation is further structured by separating quantum and 
thermal contributions as well as contributions of the low- 
est and of higher Landau levels. 

1. Low temperature approximation and separation into 
thermal and quantum contributions 

For t <C 1, we can use the asymptotic expansion for 
the digamma function and approximate the inverse fluc- 
tuation propagator, Eq. (11), as 
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Class of diagrams 


terms in Kei. 12, ci. Jiq. (Adj 


after integ. by parts 


terms in Usadel technique, cf. Eq. (Al) 


<^ 


c {sh,2c) _ . 
""DOS ^ 


r (Sh,2c) , 

""DOS 












1/1 lxAsh,2c) 
2^ 2"'^ DOS 


r (1) r {2a) 

—dale — oalc 


Sffsc 


< 


w\ 




r (cth) 


/fX (1) 




.,r (2..) f2f,! 




3B 


— ZOCTsc 


A 


= 


<x> 


-A 


-A 




MT 


Sffan 



Table I. A comparison of the expressions for the fluctuation conductivity in a jiarallcl magnetic field obtained from a diagram- 
matic calculation in Ref. 12 and from the Usadel equation approach in this manuscript. The LHS shows the classes of diagrams 
evaluated in Ref. 12 together with the corresponding terms given in (A3). The RHS shows the expressions derived from the 
Usadel equation technique as given in (Al). Aligned adjacent cells or blocks of cells are (in total) equal. 



Considering further Eqs. (3) to (9), we notice 

that terms with an integrand containing B' (uj) = 
i9<^ coth (cti/2T) only contribute at finite temperatures as 
the integrand vanishes exponentially as the temperature 
goes to zero. The terms containing B, found from 6(7dos 
and 6aic\ contribute at finite temperatures as well as 
in the zero-temperature limit. For these terms, a useful 
separation can be achieved by writing 

S(w) = [^(w)-sgn(w)]+sgn(w), (78) 

where the combination in square brackets vanishes expo- 
nentially with temperature. The notion quantum terms 
will be used for those contributions originating from 
sgn (w) in this decomposition and they will be denoted as 
Sao; they are temperature- independent. The remaining 
terms, which encode the temperature-dependence and 
originate either from B' or from the difference B (w) — 
sgn (w), will be referred to as thermal terms and denoted 
as dar; they vanish in the zero-temperature limit. 



equations A 1 and A <i; 1 define the regions I and II, 
respectively, as introduced in Sec. TIB. 

If A is large, the lowest Landau level is the only strongly 
contributing mode. The divergence close to the critical 
line is accounted for by the term n = in the sum, a re- 
striction to only this term is known as the lowest Landau 
level approximation. This approximation has been put 
forward by Galitski and Larkin^^ for the perpendicular 
magnetic field case z = 0, and allowed them to obtain 
closed formulas. 

If A is small, the significance of the lowest Landau level 
is lost. To correctly understand the cross-over, we sepa- 
rate each contribution in region I into four parts: First, 
each sum is split into the lowest Landau level (n = 0) 
term (LL), and the sum over higher Landau levels (n > 0, 
HL). Furthermore, it is of calculational advantage to per- 
form a separation of each of the resulting terms into the 
thermal (T) and quantum (0) contributions as explained 
above. 



2. Summation over Landau levels 

The evaluation of the formulas for the fluctuation cor- 
rections requires a summation over the Landau level in- 
dex n which in general is difficult to treat. However, we 
note that when approaching the transition, Lq = 1/£q 
diverges for a; = while the other L„ remain finite. In 
order to quantify the importance of the n = term, we 
introduce the parameter 

Li ln(l + /i)' 

The number A measures the relative importance of the 
lowest Landau level as compared to higher ones. The 



A 



(79) 



C. Evaluation of integrals in region I 

1. Thermal terms 

Lowest Landau level (LL): For the evaluation of the 
thermal part, only frequencies lo c:^ T <^ Tco contribute, 
and thus (77) can be expanded in lo and h. In order to 
find all relevant contributions, for the dos and an contri- 
butions it is necessary to expand to second order. The 
result of this expansion can be written as 



i5cr' 



'T,LL — ^ 



duj 



a,(B-l)uj- l3,B'uj^ 
{2aohf +LU^ 



(80) 
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where i enumerates different contributions. The prefac- 
tors ai and /3j are z-depcndcnt: 



^dos — 
(2o) 

ale = 

a^'') - 



1 

' 1+z ' 



/3. 



1 

3+2 ' 

0, 
0, 



dos 



/3(2b) _ 
h'sc — 



1_ 

1+z ' 
2 

1+z' 

0, 

1 

3+z' 
4+2z 
3+z • 



(81) 



The integrals can be evaluated as a sum over poles of B 
and B'. The result is the expression for 6aT,LL presented 
in Eq. (30). 

As compared to the result of Galitski and Larkin,^^ 
we have modified by subtracting ln(/i), which is part 
of the quantum component [it will be treated separately 
below] . and we added the angular dependence via z. Care 
must also be taken as Galitski and Larkin employed the 
traditional classification based on diagrams, while we use 
here the classification introduced in Rcf. 20. 

Higher landau levels (HL): The quantum contribu- 
tion due to the higher Landau levels are not singular 
at the transition and can be neglected. When going to 
smaller angles, the thermal contribution due to higher 
Landau levels increases. The point where the higher lev- 
els start to play a role marks the onset of the cross-over 
between region I and region IL 



2. Quantum terms 

Lowest Landau level (LL): The quantum part of the 
lowest Landau level contribution consists of terms origi- 
nating from 5a dos and 5a^^} . Contributions to the inte- 
gral are not restricted to lo < T now, thus linearization 
of the integrand is not appropriate here as it is for the 
thermal part. In these terms, w only appears in £, thus 
we can make a change of variables iuj — > Cj. The integra- 
tion contour of Cj can then be rotated by 90 degrees onto 
the real axis, thus making the integrand real. 

Turning first to the quantum part of the DOS correc- 
tion (4), it can be written as 



0,LL 



e 

2^ 



l + h 
l + z 



^ dx 

111 X 



(82) 



and the integral can be expressed in terms of the loga- 
rithmic integral function^^ Li (x) = J^^ resulting in 

the first term in Eq. (34). The second contribution to 

the quantum part originates from aic^ of Eq. (7). Upon 



setting a 



2(l+/i) 
1+z 



it can be written as 





'1 1 




/ dx 






ll+h 


X X + a 





1 

In a; 



In a; 



This gives the second term in Eq. (34). In order to ap- 
proximate this integral, we note that in the limit h <^ a, 



the asymptotic behavior at the phase boundary is due to 
the singularity of the first inverse logarithm: 



Scroti = 



27r2 



1 



Inh (/i<C a) 



(84) 



Furthermore, the next term in this asymptotic expansion 
around small h is an a-dependent constant. For more 
general values of h, the integral needs to be evaluated 
numerically. 

Higher landau levels (HL): The quantum contribu- 
tions from n > are non-singular at the transition. Still, 
one needs to be more careful than for the thermal contri- 
bution, because the sum over all levels is in fact divergent, 
and a cut-off has to be introduced. We convert the sum 
over n in the expressions (4) and (7) into an integral. 
This transformation becomes exact in the limit ^ ^ 0, 
and it is a good approximation otherwise. 

The resulting integrals are doubly-logarithmically di- 
vergent, so a cut-off A must be introduced, as discussed 
before Eq. (39). After introducing dimensionless inte- 
gration variables and again rotating the contour for w, 

we arrive at the integrals Sa^f'^l and SaQ^'^^^ given in 
Eq. (39). 

These integrals can also be expressed explicitely in 
terms of the logarithmic integral Li(a:). The correction 
then takes the form 



'0,HL = 



e 

2^ 



KLi 



1 



■ln^^+TZ(a^K^)-2n{a,K) 
In a ^ ' 



where we introduced the function 

1\ /I 



n{x,y) = X 



Li 



Li 



(85) 



(86) 



and abbreviated a = (1 + h)"^ and K = A/(2aco), and A 
is the energy cut-off up to which superconducting fluctu- 
ations are taken into account. As discussed in Sec. II B, 
it is of the order of the elastic scattering rate of electrons, 
1/r. 



D. Evaluation of integrals in region II 

For region II, the lowest Landau level loses its special 
significance. The summation over Landau levels can be 
replaced by an integration. For the quantum part Sao^HL, 
the result was already given in the previous chapter in 
Eq. (85). We thus only need to consider the thermal term 
5(Jt,hl here. 

In contrast to region I. one needs to integrate over Lan- 
dau levels also for the thermal contribution, an expansion 
of f in w and h, however, is again possible. After apply- 
ing this expansion to all thermal terms, one finds that 
the term (Jcric'''' is dominant over all other terms, which 
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axe reduced by the factor T/Tcq, which in the considered 
regime is small. This is in contrast to the case of the LL, 
where all terms contribute with equal magnitude. 

The dominant term (5(Tsc'''' can then be written as 
^F{r}), with F{ri) given by an integral over q and w/^ 



J (2„) sinh !/J(, + t2)!+j: 



2 • 



(87) 



Here, r] is given by 77 = (q;o(/i) ~ (^c)/T, which for low 
temperatures can be written as r? = 5^ \2h + /i^] . 

The integration in y can be performed analytically af- 
ter writing 1 / sinh^ y as a sum of its poles, resulting in 
the following expression 



1 f°° 

^-^ \ dx (ttx- 



V) 



i^'(x)+<(x) + ^ 



The integral F can further be simplified to give Eq. (41). 

The asymptotic behavior of F for large and small values 
of rj as stated below Eq. (42) can be found by inserting 
the asymptotic expansions i^'ix) « for small x and 
ip'ix) = x~^ + x~^/2 + x~^/6 H for large x. 



VI. CONCLUSION 

We studied the fluctuation conductivity of disordered 
superconductors subject to a magnetic field in the metal- 
lic phase using a quasiclassical kinetic equation approach. 
The derived expressions generalize the results of Ref. 20, 
in which calculations were performed for films in per- 
pendicular magnetic fields, to a more general class of 
pair-breaking transitions in superconductors of different 
geometries. We were also able to make contact with pre- 
viously derived formulas for films and wires in parallel 
fields. For the parallel magnetic field case we performed 
a detailed comparison between the traditional classifica- 
tion based on diagrams with the classification based on 
the Usadel equation approach, see Tab. I. As an ap- 
plication, we studied fluctuation corrections in a super- 
conducting film subject to a tilted field with emphasis 
on the low-temperature regime where these films display 
the phenomenon of the NM. The calculations performed 
in Ref. 20 and in this manuscript clearly show in which 
way different physical mechanisms contribute to this phe- 
nomenon: The reduced quasiparticle density of states 
leads to an increase in resistance, while Cooper pairs do 
not efficiently transport charge in this regime and are 
th^is unable to compensate this effect. 

It should be noted that besides the fiuctuation correc- 
tions additional quantum corrections exist in disordered 
electronic systems. Both the weak localization correction 
and the Altshuler-Aronov interaction correction are log- 
arithmically divergent at low temperatures. These cor- 
rections do not become singular near the phase transi- 
tion, nevertheless they need to be taken into account in 



a quantitative comparison to experimental data. The de- 
phasing time has been introduced phenonienologically 
in this paper, see Sec. II A. For vanishing magnetic fields, 
dephasing is necessary in order to regularize the anoma- 
lous Maki-Thompson correction. At low temperatures 
and finite magnetic fields, no regularization is required, 
but dephasing may still infiuence the magnetic field and 
temperature dependence of the fluctuation corrections. 
For example, it leads to a change in the phase boundary. 

The phenomenon of the NM is of particular rel- 
evance in the context of the magnetic field-tuned 
superconductor-insulator transition observed in thin dis- 
ordered superconducting films. '''^ In these films, one finds 
a change from a superconducting to an insulating behav- 
ior {dp/dT < 0) as a function of the magnetic field which 
leads to a very pronounced resistance maximum at low 
temperatures. For suitably prepared films, the resistance 
at the maximum may exceed the normal state resistance 
by many orders of magnitude. Unfortunately, transport 
in the highly resistive phase is very difficult to describe 
theoretically. The NM predicted^ ^ and also observed^'^^ 
in low-resistive films may be viewed as a precursor of 
this effect. Importantly, for low-resistive samples, con- 
trolled calculations can be performed, as was done in this 
manuscript. One may hope that a comparison to experi- 
mental data may contribute to a better understanding of 
the phenomenon of the NM and of the properties of the 
thin films in general. 
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Appendix A: Comparison to diagrammatic 
perturbation theory 

In this appendix wc show how the results of the dia- 
grammatic calculation of Ref. 12 and those obtained in 
the Usadel equation approach can be transformed into 
each other. The systems under study are films or wires 
subject to a parallel magnetic field. 



1. Fluctuation corrections obtained from the 
Usadel equation technique 

In Ref. 12, the authors studied the fluctuation conduc- 
tivity of systems for which the fluctuation spectrum is 
continuous in one or two dimensions and strongly dis- 
cretized in the transversal direction(s). For these sys- 
tems, the eigenvalues of Eq. (1) are given by Eq. (15). 
The authors assume that only the lowest mode, n = is 
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of relevance. The corresponding eigenfunctions are plane 
waves in the unconfined directions, tpq = e^'^'", where q is 
a one- or two-dimensional wave vector. 

Under these assumptions, one can also obtain expres- 
sions for the fluctuation corrections from the general re- 
sults presented in Sec. II A. Inserting the matrix ele- 
ments of Eq. (12) and Eq. (13), restricting the sum over 
transverse modes to n = and writing a±o = a, the 
expressions for the fluctuation corrections can be written 
as 



2De^ I ^^^B'Re£'lmL-Blm{£"L), 



27r (27r) 

dco d'^p B' 
27r (27r)^IV72 + Q!' 

du! d'^q 



27r (27r) 
2 2 f duj d'^q 2 a' 



2tt 



(2^) 



^qiBRe{£"£'L'), 
qlB'ReS'Re{L'^£'), 



Saif^ = 16D^e^ I ^^qlB'Im£'Re{£'L)lmL, 



27r (27r)" 



where q is the d-dimensional wave vector along the un- 
confined direction(s) {d = 1 for wires, d = 2 for films), 
we remind that B{u>) = coth(a;/2T) is the bosonic distri- 
bution function. 



£q{0j) 

: In — + 



1 Dq^ + 2a + iuj 

2 ^ 4^ 



,(A2) 



and L = l/£ is the fluctuation propagator. In the dis- 
played formulas the momentum and frequency arguments 
{q, ui) have been suppressed for the sake of brevity. 



2. Fluctuation corrections obtained in Ref. 12 

We next compare the results (Al) to the results derived 
diagrammatically in Ref. 12. In that publication, the 
the following fluctuation corrections are presented [we 
adjusted them to the notation used in this manuscript] 



j: (sh,2c) 



c (sh,3c) 



'AL 

(cth) 



' AL 



2 f d'^'qdcdu) , 



{2ttY2-k 



n (e) B' {to) L; {to) Re [DqlC^^{2e - u;)] 



= AiD'e' f (e) B'{^)L;{c.)C, (2e - uj) C, (-2e + a;) , 



2A = -16De' [ ^;^H (e) B {co) [0^ (2e - u;) L* (c.)] , 

J (27r) ZTT 



B = -72De' [ (e) B (u;) L; (c.) DqlC'^ (2e - c.) , 

J (27r) 27r 



Sai''^^ = -8De 



(27r) 
2 f d'^qduj , 



B'{uj)Dql (Re[L,(a;)5;(a;)])^-Im L,{uj) {£'^{uj)y ImL,(a;) 



(27r)''27r 

/ (S^^ ^'^^ [L^(c.)f;(a;)f^(u; 



(A3) 



The corresponding diagrams are shown in Tab. I. 



that 



M'J=Sa^fj^'\ (A4) 



AL 



(A5) 



The first five terms in Eq. (A3) are not yet integrated 
with respect to the fermionic frequency e, while the 

last two terms are already very similar in structure to The integration in e for the first five terms in Eq. (A3) can 
daii\ dai'l"'^ and Sai'l''^ in (Al). In fact, one can verify be performed analytically. The required manipulations 
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are similar to those described in App. B. The results are 
^<6T = 4e^i? / I^B' Imi ReS' , (A6) 



-2A = -8e^Dl^Blm[Ln, 

9 
2 



lB = -12e^Djj^BD4S.lLe-]. 



Comparing these expressions to (Al), we find: 

2Sados = Sa^i'liT - A, (A7) 

(A8) 



Finally, by transforming the momentum integrals to 
spherical coordinates and integrating by parts in the ra- 
dial variable, one can show that 



(A9) 
(AlO) 



2""sc — ""DOS ""DOS ' 

0=lA-3Sai'J -Ib. 
2 ""2 

We note that for dimensions d > 2 some of the terms 
in Eqs. (3) to (9) require an ultraviolet regularization. 
Indeed, an upper cut-off \Dq^ < A for the inoniciitum 
integral needs to be introduced as discussed in connection 
with Eq. (39). In 2d this divergence is very weak (doubly 
logarithmic), in 3d it is more severe. This is not unusual 
as we work with a low-energy theory that ceases to be 
accurate at higher energies. The important point is that 
the theory captures correctly the sensitivity to tempera- 
ture and magnetic fields. It should be remarked in this 
context that the boundary term, which we dropped when 
performing the integration by parts to obtain Eqs. (A9) 
and (AlO), is in fact of order A'*~^/lnA. The boundary 
term that arises once a finite cut-off is introduced, how- 
ever, is insensitive to small changes in the parameters, T 
and h. We therefore do not attribute particular impor- 
tance to the difference in the Tiltraviolet regularization of 
the terms obtained from the Usadel equation and from 
the diagrammatic technique. 

The identities (A4), (A5), (A7), (A8), (A9) and (AlO) 
show that our results, Eqs. (3) to (9), when applied to 
the case under study, are equal to the corrections derived 
diagrammatically, Eq. (A3), in one spatial dimension. In 
2d and 3d, the equivalence still holds up to details of 
the ultraviolet regularization. The identification of cor- 
responding terms is summarized in Tab. I. 



Appendix B: Frequency integration 

The derivation of Eqs. (3) to (9) from Eqs. (66), 
(67) and (68) as well as the derivation of Eq. (A6) 



from Eq. (A3) both involve an integration over the 
"fermionic" frequency e. By way of example, we pro- 
vide some details on the evaluation of the integrals. The 
integrand typically comprises two Cooperon propagators 
C„(e) = (ie — 2a„)~^, multiplied by one or more fac- 
tors of the fermionic equilibrium distribution function 
He = tanh(e/2T) or its derivatives. 

Integrals with two Cooperons can be treated by a par- 
tial fraction decomposition. 

Anm = J^-HeCn (2e - w) C„ (2e - w) 

= \ X / ne[Cr,{2e-uj)-Cm{2e-co)] 

5* (lu) - (oj) 



Aiir {an - am) 



(Bl) 



where we used the shorthand notation f = f 4^. The 
integral has been performed by the method of residues, 
the poles of "He lie on the imaginary axis at e„ = 
AniT {n+1/2) , with n integer. The resulting sum can 
be expressed in terms of the digamma function, and fur- 
ther in terms of £, Eq. (11). For am = an, one may use 
the relation da„£n ('^) = —2idu;£n (w) and finds 



1 

'2n 



CM- 



(B2) 



Integrals involving He-u can be reduced to the dis- 
cussed examples by a shift in the integration variable 
and complex conjugation, e.g.. 



7^,_.C2(2e-a;) = -^:„. (B3) 



Another useful relation, that can be obtained in a similar 
way, is 



He Cni2e-uj)Cmi-2e + uj) 



£*n (^) - £ra (^) 

Ai-K {an + am) 



(B4) 



Finally, integrals including H'^ can by reduced to the for- 
mer integrals using integration by parts, yielding 

Bnm = ^ K Cn (2e - uj) Cm (2e -uj) = 2A;„, (B5) 
Bnm = ^ K Cn (2e - uj) Cm (-2e + uj) = 2A'„^. (B6) 

If the integrand contains more than one factor of T-Le or 
He-, the following two identities can be used to simplify 
the expression 

UeV-e-^ = l-B^{He~ ^e-c) , (B7) 
H'e'He-u} = —Buj'H-e ~ ^'u (^e ~ "H-e-ui) ■ (B8) 

The first identity is essentially the addition theorem for 
tanh'^^, while second identity is obtained from the first 
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one in three steps: A shift e e + w is followed by is a step required for obtaining Eq. (4) from Eq. (66) 
a differentiation with respect to uj and a second shift 

With the stated integration formulas at hand, in com- Je 
bination with relations (B7) and (B8), all the required in- = —BuBnn — ^'^i-^nn + -^nn)- (-^9) 

tegrals can be performed. Let us give an example, which 
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